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ABSTRACT 

Context. 
>0 , Aims. Interstellar dust grains, because of their catalytic properties, are crucial to the formation of 

H2, the most abundant molecule in the Universe. The formation of molecular hydrogen strongly 
depends on the ability of H atoms to stick on dust grains. In this study we determine the sticking 
coefficient of H atoms chemisorbed on graphitic surfaces, and estimate its impact on the forma- 
tion of H2. 
CL(' Methods. The sticking probability of H atoms chemisorbed onto graphitic surfaces is obtained 

Q . using a mixed classical-quantum dynamics method. In this, the H atom is treated quantum- 

■j-J , mechanically and the vibrational modes of the surface are treated classically. The implications 

C/2 . 
C^ , of sticking for the formation of H2 are addressed by using Kinetic Monte Carlo simulations that 

follow how atoms stick, move and associate with each other on dust surfaces of different temper- 
ature. 
> 
fT^ ' Results. In our model, molecular hydrogen forms very efficiently for dust temperatures lower 

O^ ' than 15 K through the involvement of physisorbed H atoms. At dust temperatures higher than 

1/^ ' 15 K and gas temperatures lower than 2000 K, H2 formation differs strongly if the H atoms com- 

r — ^ . ing from the gas phase have to cross a square barrier (usually considered in previous studies) or 

a barrier obtained by DFT calculations to become chemisorbed. The product of sticking times 

efficiency can be increased by many orders of magnitude when realistic bamers are considered. 

If graphite phonons are taken into account in the dynamics calculations, then H atoms stick better 
(yj " on the surface at high energies, but the overall H2 formation efficiency is only slightly affected. 

$_( ' Our results suggest that H2 formation can proceed efficiently in photon dominated regions. X-ray 

dominated regions, hot cores and in the early Universe when the first dust is available. 

Conclusions. 

Key words. ISM: dust, extinction - ISM: molecules - ISM: abundances 

1. Introduction 

Interstellar dust grains play a very important role in the chemistry of the interstellar medium 
(ISM). Because of inefficient gas ph ase paths to form H2, d ust grains are considered as the fa- 



o 



vored habitat to form H2 molecules (IGould & Salpeterlll963l) . In the past decades, a plethora of 



laboratory experiments and theoretical models have been developed in order to understand how 
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the most abundant molecule of the Universe forms. The sticking of H atoms on surfaces has re- 
ceived considerable attention since this mechanism governs the formation of H2, but also other 
molecules that contain H atoms. The s ticking of H atoms o n dust grains can also be an impor- 
tant mechanism to cool interstellar gas dSpaans & SiUaEOOOl) . As H atoms arrive on the dust, they 
can be weakly bound (physisorbed) or strongly bound (chemisorbed) to the surface. The stick- 
ing of H in physisorption state has been highlighted by several experiments on different type of 



surfaces (IPirronello et al 



1997b . 



1999, 



a 



2000llPerrv & Pricdl2003h . Experiments performed on sur 



faces at higher temperatures revealed how H atoms stick in chemisorbed state ( iZecho et al. 



Hornekaer et al 



2006, 



I 



2002, 



Mennellall2006^ . In addition, theoret ical studies confirmed that the st i cking 



of H a toms varies strongly with gas and dust temperatures JLeitch-Devlin & Williamslll985 



1989, 



Shaet 



Morisset et al. 



Medina & Jackson 



2008, 



Morisset & Allouche 



2008, 



Cuppen et al. 



Buch 



201c . 



20051 
20101) . To study the sticking of an atom on a surface, it is necessary to take into 



account the vibrational modes (phonons) of the surface (iBurke & Hollenbach 
Morisset & Allouchell2008 , 



1983 



Buch 



1989 



Morisset et al .1120 101 ) 



In the ISM, dust grains are mainly carbonaceous particles or silicates, with various sizes, and 
a large fraction of the available surface for chemistry is in the form of very small grains or pol- 



yaromatic hydrocarbons (PAHsi lWeingartner & Drair 



i 



graphite surfaces , as shown by different c alculations dJeloaica & SidisI 



2001 ). PAHs have similar characteristics as 



1999, 



FerroetalJ2008) 



whereas 



Ferro et alj ( 120081) and 



Sha&Jacksonll2002 , 



3E 



Jeloaica & Sidis ( 1999 ) studied the i nteraction between H atoms and a coronene. 



Sha & JacksonI ( 120021) determined the H-graphite interaction. All 



these theoretical works obtained identical characteristics: the H atom coming from the gas phase 
can physisorb (43 me V) or chemisorb (0.76 eV) on the gr aphite or PAH. Several experimental stud- 



ies sh owed that H atoms can physisorb on c arbonaceous ([Pirronello et al.lll997bl l 

cates ( IPirronello et al. 1997a) and graphi tic ( Perrv & Pricell2003b surfaces, and can also chemisorb 



on graphite surfaces (IZecho et al. 



2002, 



Hornekaer et al. 



1999 . 



20001), sili- 



20061) . These experiments revealed that 



adsorption of H atoms (physisorption and chemisorption) is of key importance to form H2 for a 
wide range of dust and gas temperatures, as observed in the ISM. The H atom can chemisorb only 
on top of a C atom and a surface reconstruction is observed: the C atom goes out of the surface 
plane. This phenomenon is called puckering, and creates a barrier against chemisorption of .2 eV 



i 



I 



dJeloaica & SidisI 



1999, 



Sha&Jacksonll2002, 



Ferro & Allouchell2003 



Hornekaer etal.112006). This 



barrier is important and does not allow an efficient sticking mechanism for atoms coming in at low 
energies. Therefore, sticking through chemisorption strongly depends on the energy of the incom- 
ing atoms. Also, the sticking is sensitive to the temperature of the dust since the surface excitation 
modes play a role in redistributing the excess energy of an atom. 

The H2 formation occurs through Eley-Rideal (ER) or Langmuir-Hinshelwood (LH) mech- 
anisms. In the ER mechanism, an H atom coming from the gas phase collides with an H atom 



which is initially adso rbed (chemisorbed) on the graphite surfa ce. Theoretical dRougeau et al 



2006, 



Ferro et al. 



20081) and experimental dHornekaer et al 



2006) studies show that H atoms are 



chemisorbed on the surface in dimer configurations, and that two hydrogen dimers configurations 
are stable: the ortho-dimer and the para dimer (Figure 3, right panel). Density Functional Theory 
(DFT) calculations show that while the chemisorption of one H atom is associated with an 
important barrier, the formation of the para dimer is barr ier-less and that the formation of the 
ortho dimer has a reduced barrier dRougeau et al.l d2006h ). The formation of H2 that involves 
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chemiso rbed atoms through the ER mechanism has been studied by Den s ity Fu nctional Theory 



(DFT) (iJeloaica & Sidis 



2004b . 



1999 



calculations (Rutigliano et al. 



Sha & Jackson 



2001, 



Ree et al. 



2002 



200: 



J 



Ferro & Allouche 



Sha & Jackson 



2002, 



20031) and dynamics 



Morisset et al. 



2003 . 



Martinazzo & Tantardinill2006al) . Based on these calculations, different mechanisms have 



been proposed to contribute to the H2 for mation through the ER mechanisms: the direct ER 
that involves isolated H atoms (monomers, 



Sha & Jackson 



2002, 



Morisset et al. 



2003 



2004b, 



Martinazzo & Tantardini 



states (iBonfanti et al 



2006a), barrier-les s formation of H2 involving one H atom in a para- 



dimer configuration ( Bachellerie et alj|2007h and formation by diffusing H atoms in physisorbed 



20071). In the LH mechanism, the two H atoms are adsorbed (physisorbed) 



on the graphite surface, diffuse on the surface and collide to desorb in a hydrogen molecule. Thi s 



mechanism has also been studied (IMorisset et al 



2004a . 



2005, 



Martinazzo & Tantardini 



2006bh . 



The different mechanisms to form H2 on graphite surfaces are the following: 

Ml) LH mechanism: Two physisorbed H atoms encounter each other on the surface; 

M2) direct Eley-Rideal mechanism involving H atom chemisorbed in a monomer (only one H 

atom chemisorbed on the cycle, fig. [3] right panel); 

M3) direct Eley-Rideal mechanism involving an H atom chemisorbed in a dimer (H atom 

chemisorbed in dimer position on the first cycle, fig.[3j right panel); 

M4) fast diffusion of physisorbed H atoms that enter chemisorbed sites occupied by H atoms 

(monomer); 

M5) fast diffusion of physisorbed H atoms that enter chemisorbed sites occupied by H atoms 

(dimer); 

M6) direct Eley-Rid eal mechanism involving an H atom in meta, 5 and S (second cycle) positions 



(iDumont et al. 



20081 see figure 3 right panel); 



M7) ER mechanism by fast diffusing H atoms in the physisorption state ( IBonfanti et al.ll2007h . 



Recent experimental studies dLemaire et al.l l2010i llslam et al.l |2007|) show that the ro- 
vibrationnal states of the H2molecule formed on the surface can be detected (siUcates and graphite). 
The ro-vibrationnal distribution of the H2 allows to classify the LH and the ER mechanisms to form 
the molecule as function of the dust temperature. Experimentally, it is not possible to identify the 
Ml to M7 mechanisms proposed to form H2. The objective of our work is to classify the pre- 
dominant mechanisms to form H2 as a function of gas and dust temperature. In each mechanism 
proposed, the sticking of the H atoms in physisorption sites and chemisorption sites (figure 3 left 
panel) is crucial to allow the H2 formation. In this study, we concentrate on understanding how the 
first H atom chemisorbs on the grain with a direct ER mechanism and how the sticking probability 
impacts the H2 formation. During the collision between an atom and a surface, the collisional en- 
ergy is distributed between the vibrational excitation of the newly formed bond and the vibrational 
modes of the surface. Therefore, to calculate the sticking of atoms on a surface in a realistic way, 
one has to consider the lattice dynamics. For this purpose, we study the interaction of H atoms with 
graphitic surfaces, which are surfaces having similar properties as PAH surfaces. Our attention fo- 
cusses particularly on these su rfaces because PAHs represe nt a large fraction of the total surface 
area of dust grains (up to 50%. IWeingartner & DrainduOOlb . 

In the first section, the time dependent dynamics method is presented to study the sticking of 
H atoms in chemisorbed sites on the graphite surface in the collinear geometry. In the subsequent 
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section, the different mechanisms for the formation of H2 are implemented into the Kinetic Monte 
Carlo (KMC) simulations. The sticking probability obtained by the dynamics method is included 
in the KMC code. In the last section, we assess how this sticking changes the processes involved 
in the formation of H2 and identify the predominant mechanisms to form H2 molecules as function 
of the gas and dust temperature. 

2. Sticking: Dynamical calculations 

In this section, the time dependent wave packets propagation method is presented to calculate the 
transmission probability for an H atom to overcome the chemisorption barrier The time dependent 
dynamics method, that takes into account vibrational modes of the surface, is used to calculate 
the sticking probability of H atoms in chemisorbed sites on a graphite surface. Our calculation are 
performed in the collinear approach. 

2.1. The wave packet propagation method 

This method solves for the time dependent Schrodinger equation: 



d \Mt)) 
ot 



(1) 



where H is the Hamiltonian of the H-graphite system which is the sum of a kinetic operator and a 
potential operator and |i/'(0) is the time-dependent wave function. This method consists of choosing 
an initial wave function (at time t=0) to model the H atom in the gas phase. The application of the 
evolution operator on the wave function allows us to obtai n the wave function at time t + At. 



The time propagation is performed using a Lanczos method (iLeforestier et alJ 



1991). The kinetic 



19831) . The 



operator is applied on the wave function using a Fourier meth od jKosloff & Kosloftl 
propagated wave function is analyzed by a flux analysis method (iBalint Kurti et alJll990h to obtain 



physical information such as the probability of transmission of a potential barrier. 



2.2. The probability of transmission of the chemisorption barrier 

The interaction potential between the H atom and the graph ite surface is well known 



(lFerro&Allouchd2003 . 



Sha& Jackson 2002. 



Jeloaica & Sidis 



19991) . This potential is represented 



on the figure [T] as a function of the reaction coordinate z between the H atom and the graphite sur- 
face in the collinear approach. The main characteristics of the p otential derived from t he DPT 
code PWSCF included in the QUANTUM ESPRESSO package of Joiannozzi et al.i2009L are flie 
foUowing: 

1 . A surface reconstruction is observed: the H can be chemisorbed on the graphite surface on top 
of a C atom with an energy bond of 0.76eV. The C atom where H is chemisorbed, goes out of 
the surface plane. 

2. Chemisorption on the surface has an associated barrier of 0.2 eV. 

The probability of transmission (FiglT]i of the chemisorption barrier is calculated using the 
wavepackets propagation method (described in section l2Tt as a function of the incident energy of 
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the H atom coming in from the gas phase. The probability of transmission decreases with decreas- 
ing incident energy of the H atom. For energies smaller than 0.2 e V (the value of the chemisorption 
barrier), the H atom gets through the barrier by the tunneUng effect. 

The DFT calculations are performed on a 2x2 working cel l. In this case , the ch emisorption well 
is about 0.76eV and the chemisorption barrier is about 0.2eV.|( 



Casolo et al. 



(2009) have shown that 



the chemisorption well is affected by the size of the working cell. For example, for a 5x5 working 
cell, the chemisorption well is 12% higher than with a 2x2 working cell and the chemisorption 
barrier is of 0. 15 eV. In our dynamics calculations, we consider a chemisorption barrier 0.2eV. The 
value of this barrier will affe ct the sticking prob abilities and the transmission probability of the 



barrier For a smaller barrier dCasolo et al. 



20091) the sticking probabilities and the transmission 



barrier calculated using dynamics method would be higher than the values calculated here in the 
energy domain between 0. 15eV and 0.2eV. 



2.3. Sticking of H on grapiiite using a mixed ciassical-quantum dynamics metliod 

The detail of the theory to study the sticking of H atom on a graphite surface using the mixed 



Morisset et al 



(I2OIOI) . Here, we 



classical-quantum dynamics method is presented in the paper of 1] 

present a brief description. The hamiltonian used to describe the interaction of the atom with the 

dynamical surface is: 



H^T + Vtot + H^, 



(2) 



where T is the kinetic operator, Vjot i^ the gas-surface potential operator, H^ is the Hamiltonian 
for the lattice vibration without coupling to an atom. In the mixed classical-quantum method, it has 
the following form: 



'■=' j=l 



(3) 



where A^ and M are the number of phonon modes and phonon bands, respectively. Q. . and P. . 
are the positions and momenta of the vibrational modes of frequency wj; . In the previous papers 
of Morisset & Allouche (2008, 2010), the Taylor series expansion of the potential V,oi in terms 
of the displacements of the lattice atoms allows one to write the potential as the sum of the static 
potential Vo{z) between H and the graphite surface in the collinear approach plus an interaction 
potential between the atom and the bath of phonons. The Taylor expansion is truncated to the 
linear term in phonon coordinates, which is the one-phonon exchange approximation. In the mixed 
classical-quantum approach, the total gas-surface interaction can be written as: 

N M 

V,o, = Voiz) + YjYj SiA^^QiJ' (4) 

'=' j-1 

where gij{z} represents the coupling terms between the motion of H and the bath of phonons 
with a vibrational mode of frequencies w/y and a polarization vector eg. To calculate from DFT 



Morisset & Allouchd (l2008b . The 



calculations the gij{z) terms, a model has been developed by 

giiiz) terms are directly ca lculated of DFT phonon calculation in the harmonic approximation 



dAshcroft & Mermin 



19761) . These calculations allow to obtain the dynamical matrixes D in the 



irreducible Brillouin Zone (IBZ). By diagonalization of D, the frequencies and polarization vector 
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of the vibrational mode are ob tained. The phonon di s persio n of the H-graphite system is presented 
in the Figure 2 of the paper of lMorisset & Allouchd (l2008h along the high- symmetry directions in 
the IBZ. The phonon dispersi on obtained by DF T phonon calculations are in good agreement with 
previous paper on the subject dMohr et alj|2007l) 

The dynamics method is performed through a mixed classical-quantum approach. In this 
method, the H motion (z coordinate between H and the graphite surface) is treated quantum- 
mechanically using the wavepackets propagation method (described in section l2Tb . whereas each 
vibrational mode is treated classically. The time dependence of the classical variables is obtained 
by solving the equations of motion derived from the Hamilton equations; 



Ql 



dH, 



eff 



PiJ = - 



dPij ' 
dH.ff 



dQu ' 



(5) 



where H^ff is the effective Hamiltonian, given by the expectation value of the Hamiltonian H : 



H, 



eff 



--<^\H\il,> 



(6) 



Initially, each vibrational mode obeys the equipartition theorem jMcOuarridbOOOl) : there is -^ of 
energy in each vibrational mode where T is the surface t emper ature. The wave function is analyzed 
using the reactive flux as described in 



Morisset et al 



(120 101) . which allows one to calculate the 



sticking probability as a function of the incident energy of the H atom and the surface temperature 
T. 



2.4. Dynamics results 

2.4.1. Transmission coefficients 

In previous studies, we calculated the probability for an atom to stick on a grain surface to be 



equal to the transmission coef ficient to cross the barrier against chemisorption JCazaux & Tielens 



2002, 



a 



Cazaux & Tielensll2004l) . We considered the barrier to be square, and obtained transmission 



coefficients that strongly depend on the energy of the incoming H atom (see figH] right panel, 
dotted lines). The transmission coefficients are very small at energies below 0.2 eV (~2320K, 
which corresponds to the height of the barrier). 

However, the barrier against chemisorption is far from being square. Actually, DFT calculations 
show that this barrier varies from 0.7 9 A to 1 .05 A in width, with a height of 0.2 e V (see figlH 
left panel. iMorisset & Allouchdl2008h . The transmission coefficients obtained with such a barrier 



are shown as solid lines in figH] right panel. It is clear that the approximation of a square barrier 
underestimates the transmission coefficients by orders of magnitude. Already at energies of 1000 K, 
an H atom would have 10^ times more chance to cross a realistic barrier as calculated by DFT 
than a square barrier. This effect has profound consequences for the sticking of H atoms at low 
gas temperatures, and will influence which processes govern the formation of H2 and with what 
efficiency. 



S.Cazaux et al.: When sticking influences H2 formation 



2.4.2. Sticking 



In figlH we present the sticking of H atoms on graphite surfaces calculated with the mixed classical- 
quantum dynamics method in order to consider the dissipation of the H atom's energy through the 
diiferent modes of the lattice. In the present calculation, N=ll phonon bands with M=61 vibra- 
tional modes have been included; the phonon bands of C-H (H chemisorbed on the graphite), the 
four phonon bands of C-C of higher frequencies, the longitudinal mode of C-C, and the three acous- 
tic bands of C-C. Our dynamics calculation shows that the energy exchange between H and the 
vibrational modes of the substrate takes place by acoustic phonon modes of C-C. The chemisorp- 
tion barrier (0.2 eV=2320 K) governs the sticking mechanism: the sticking probability decreases 
with the incident energy for each surface temperature T^,,,,. At incident energies smaller than the 
chemisorption barrier, the sticking coefficient decreases with increasing Tdust- At low dust tem- 
peratures, the surface allows better dissipation of the energy of the incident atom, and therefore, 
facilitates the sticking. For example, the difference in sticking efficiency on surfaces of 10 K or 
on surfaces of 125 K can be a factor of 3 for incident energies around 1000 K. At incident en- 
ergy higher than the chemisorption barrier, the variation of the sticking probability with the surface 
temperature is negligible. 

Our calculation can be compared to experimental work by 



Zecho et al. 



(120021) . For a surface 



temperature of 150K and with a beam centered about 0.2eV, these authors obtain a sticking prob- 
ability of 0.5. In our dynamics calculation, in the same condition in temperature and energy, our 
sti cking probability is about 0.6. Our results are in good agreement with the experimental work 



of 



Sha et al 



(l2005h and 



Zecho et alj (12002 ).The other theoreti cal calculations on the subject have been performed by 



Kerwin et alJ (120061) who obtain a probability of 0. 1 and 0.0025 in the same 



domains of energy and temperatures, respectively. These probabilities are smaller because these au- 
thors included in their dynamics calculations a coordinate to simulate the energy dissipation in the 
surface, but do not explicitly take into account the vibrational modes of the substrate. Therefore, to 
study theoretically the sticking of an atom on a graphite surface, it is necessary to take into account 
the vibrational modes of the surface in the dynamics calculations. 



3. H2 formation using Kinetic IVIonte Carlo simulations 

We use a step by step Monte Carlo simulation to follow the chemistry occurring on dust grains. 
This method allows us to follow each individual H atom arriving, binding and moving on the 
surface, and forming molecules through different mechanisms. However, this method does not 
give us informations on the energy of the atoms and formed molecules. The graphite surface is 
comprised of benzenic rings that allow the H atoms to physisorb or chemisorb, as presented in 
figS H atoms can physisorb above each C atom, and on the bridge between 2 C atoms ( position E). 
Howe ver, recent studies showed that H atoms cannot physisorb at the ce nter of the r i ng ( Ferro et al 



2OO2I) . The physisorbed energy has been chosen to be 43 meV from 



Ohio et al 



Bonfanti et al. 



J 



mm and the 



20071) . H atoms 



barrier between 2 physisorbed sites is almost negligible (5-7 meV, 
can chemisorb on top of C atoms with an associated barrier of 0.2 eV and an energy of 0.76 eV. 
If an atom is already present on a ring, a second atom, with a spin opposite to the adsorb ed atom, 
can become chemisorbed in a para configuration without a barrier dRougeau et al.ll2006l) . making 
a para-dimer. H atoms can also become chemisorbed in ortho and meta configurations, but these 
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processes are associated with a barrier of 0.26 and 0.16 eV respectively (Rougeau et al. 2006)The 
para-dimer is the favorable configuration to form H2 because it is barrier-less compared to the other 
configuration. A third a tom arriving on this dimer can create an H2 molecule without a barrier 



I 



jBachellerie et alj|2007h . If an atom from the gas phase arrives on an adsorbed atom (a monomer), 



then there is a barrier of 9 meV, as derived by 



Morisset et al 



(l2004al) 



Species that are accreted from the gas phase arrive at a random time and location on the dust 
surface. This arrival time depends on the rate at which gas species collide with the grain. This rate 
of accretion can be written as: 

Race = riHVHO-S, (7) 

where Rh and V// are the density and velocity of the H atoms in the gas, cr is the cross section 
of the dust particle and S is the sticking coefficient of the H atoms with the dust. The sticking 
in chemisorbed sites has been discussed in the previous section. This coefficient strongly depends 
on the form of the barrier to chemisorb, but is also somewhat sensitive to the surface temperature 
since energy transfer to or from the phonon can be important to loose excess energy. In figure|2l left 
panel, different sticking probabilities for the physisorption case are presented as a f unction of gas 



tempe rature (Tgoj). In this figure, SI represents the sticking probability derived bv ICuppen et al 



(120 101) based on a soft cube model: 



Sphys - (1 + «! y/Tgas + T grain + a2Tgas - O^T^^^) , (8) 

with ai = 4.2 x 10"- K'^^^, aa = 2.3 x 10"^ K"', and aa = 1.3 x 10"^ K^^. SI decreases with the 
energy of the gas. At 100 K, Sl=0.8 while it becomes 0.2 at 1000 K. At low Tg„s, SI is sensitive 
to the temperature of the dust. In this model, the sticking probability of H in physiso rbed sites on 



graphite is important. S2 represents the sticking probability obtained bv lBuchI ( 119891) . This model 
is based on the quantum mechanical perturbation theory. This calculation includes a microscopic 
description of the solid structure and vibrations. S2 is very weak, and varies between 0.02 and 10 "* 
on the energy domain Tg^j^SO-SSO K. 

Atoms arriving from the gas phase on top of a carbon atom can either be physisorbed (with 
a sticking efficiency Sphys) or be chemisorbed (with a sticking efficiency Sehem)- In our model, 
because of the high barrier to access chemisorbed sites, H atoms mostly arrive from the gas phase 
in physisorbed sites for low temperature gas. 

The species that are present on the surface can go back into the gas phase though evaporation. 
The evaporation rate of H atoms in physisorbed or chemisorbed sites can be written as: 

En 

RevapiH,) = VpXeXp (--—), (9) 

Ec 

RevapiH,) = VfXeXp (--—), 

KbI 
where Ep and E^. are the binding energies of the H atom in a physisorbed or chemisorbed site, re- 
spectively, and Vp and Vp are the oscillation factors of the atoms in the physisorbed and chemisorbed 
sites taken as Vp=10'"s"' and V£.= 10'"'s"' . 

The species that arrive at a location on the surface can move randomly through tunneling ef- 
fects and thermal hopping. The diffusion rates for an atom to go from a physisorbed to another 
physisorbed site (Rpp), and for an atom to go from a physisorbed site to a chemisorbed site (Rpc) 
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can be written as: 



R 



pp 



VpXexpi—-—) 
Kb I 



"pc — Vp X rpc. 



(10) 



where Ppc is the probabiHty for the atom to go from a physisorbed site to a c hemisorbed site 
by tunneling effects or thermal hopping, as described in ICazaux & TielensI (120041) . When thermal 



hopping dominates, which is usually the case for high temperatures, or when the barrier is low, 
this probability can be written as Ppc - exp(- j-y), where Epc is the energy of the barrier between 
the physisorbed and the chemisorbed site. This barrier is 0.2 eV, which means that physisorbed H 
atoms cannot directly reach chemisorbed sites. However, if a chemisorbed atom is already present, 
then the barrier in its associated para-site disappears. In this case, the physisorbed H atom can reach 
these para-sites without a barrier and consequently will hop thermally in chemisorption. 

The mechanisms included in our KMC simulations are the mechanisms Ml to M5 described 
in the introduction. 



Dumont et al 



(120081) have included in their KMC calculations formation of 
ortho, para, meta, and 5 and S dimers (see fig [3] right panel). The H atoms can difiiise from one 
dimer position to another. These simulations explained the occurrence of two peaks in the ther- 
mal recombinative desorption of molecular H2 from clean graphite surfaces, after chemisorption 
of H atoms. The first peak (450 K) is essentially due to the para-dimer desorption whereas the 
second peak (560 K) is due to the ortho-dimer desorption. The others mechanisms included in 
their simulation contribute to the H2 desorption but with a small amount. In this work we con- 



centrate on H2 formation on dust 



grains with temperatures smaller than 150 K. The mechanisms 



considered bv lDumont et alj (l2008b are not relevant for the range of temperatures considered in this 
study since diffusion of chemisorbed H atoms as well as thermal recombinative desorption occur 
at much higher temperature s. In this study we use similar mechanisms as the ones considered in 
Cuppen & Hornekaen (120081) . and also include the formation of H2 with 2 physisorbed atoms. 

The mechanisms to form H2 are associated with very high barriers, if they involve a 
chemisorbed H atom in a monomer (mechanisms 2 and 4), and are barrier-less if chemisorbed 
H atoms in para-dimers are involved. We perform several simulations for the formation of H2 on 
graphitic surfaces with various grain and gas temperatures. We first consider the barrier to enter 
chemisorbed sites as being square, and then consider more realistic barriers obtained by DFT cal- 
culations. Finally, we also take into account the efi^ect of the phonons on the sticking of atoms in 
chemisorbed sites. 

In our simulations, we calculate the total sticking coefficient times the total H2 efficiency, Sxe. 
This product can be written: 



0X6 — pj^ys X ^phys + O clwm ^ ^chem ■> 



(11) 



where Sph^s x £phys is the product of the sticking coefficient in physisorbed sites times the H2 
efficiency involving physisorbed atoms, and Sdmn x fe/iem is the product of the sticking coefficient 
in chemisorbed sites times the H2 efficiency involving chemisorbed atoms. The efficiency of H2 
formation is written as: 

2 X hh. 



«H 



(12) 



where n//, is the number of H2 molecules formed, and Wh the number of H atoms that arrived on 
the surface. 
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3.1. Effect of the barrier on Hj formation 

The transmission probability for H atoms to overcome the chemisorption barrier (0.2 eV) in the 
collinear approach is calculated following two different approximations: 

1) the chemisorption barrier is considered as a square barrier (figure [1] right panel, dashed lines). 
The transmission probability is presented in figure[T] left panel (dashed lines), 

2) the chemisorption barrier is obtained by DFT calculations (figure [U right panel, solid lines) 
using the wavepacket propagation method described in paragraphs 2. 1 and 2.2, and the transmission 
probability is presented in figurefl] left panel (solid lines). 

The simulations for the formation of H2 considering a square barrier and a barrier obtained by 
DFT calculations to enter chemisorbed sites are shown in figlD The square barrier model has been 



used in previous studies (ICazaux & Tielensll2002l l2004l) . The objective of this study is to under- 



stand the influence of the form of the chemisorption barrier on the H2 formation using the KMC 
simulations. For the H2 formation, mechanisms Ml to M5 are included in the KMC simulations. In 
this figurelH ei represents the efliciency of H2 formation when a square barrier is considered, while 
62 represents the same efficiency when the barrier is obtained by DFT calculations. We considered 
temperatures of the gas of 100, 500, 1000 and 2000 K with grain temperatures from 8 to 90 K. 
Our results show that: At low dust temperatures (T^„jr<15 K for barrier 1 and < 20 K for barrier 
2), the formation of H2 involves physisorbed atoms, either with the mechanism Ml (2 physisorbed 
atoms) or with the mechanism M5 (one physisorbed atom entering a chemisorbed site occupied by 
an H atom in a dimer). The efficiency to form H2 at these temperatures is on the order of 100%, and 
therefore the product Sxe is equal to Siphys- At dust temperatures higher than ~15 K, the physisorbed 
atoms evaporate, and molecular hydrogen is formed with chemisorbed atoms. The type of barrier 
considered has a very important effect on the formation of H2 since it determines the sticking of H 
atoms in chemisorbed sites. At low gas temperatures, the H atoms have a low probability to enter 
the chemisorbed sites, and therefore the formation efficiency is low. H2 is formed through the di- 
rect Eley-Rideal mechanism when an H atom from the gas arrives on a dimer on the surface (M3). 
Once a first H atom sticks to one chemisorbed site, the second one sticks to a para-site without a 
barrier, and a third H atom coming on the dimer can form an H2 molecule without a barrier The 
atom that stays on the surface is used to make another dimer, again without a barrier, and again 
another H2 molecule can be formed. With this process, 5 x e is orders of magnitude higher than 
the sticking alone. As the gas temperature increases, the barrier to become chemisorbed becomes 
easier to overcome, and H2 can form through the direct Eley-Rideal mechanism when an H atom 
from the gas arrives on a single H on the surface (M2). 

3.2. Effects of the surface dynamics on the formation of Hi 

We used Kinetic Monte Carlo simulations including mechanisms Ml to M5 to calculate the for- 
mation of Ht on dust surfaces of 10 and 125 K as function of the gas temperature. The atoms 
present in the gas phase arrive on the surface with an energy that allows them either to stick in 
chemisorbed sites (if the barrier to chemisorption can be overcome), or to stick in physisorbed 
sites (if atoms can thermalize in physisorbed sites), or to bounce back to the gas phase. For the 
two surface temperatures considered, the sticking in chemisorption has been calculated previously 
(section 2.4.2), and takes into account the effect of surface dynamics. However, the sticking in ph- 
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ysisorption varies strongly from one study to another. In this work we consider strong sticking in 

J II 1 J II 1 

physisorption (ICuppen et al.1120101) . and weak sticking in physisorption jBuchll 19891) . and calculate 



for each case the product of sticking times efficiency for H2 formation, S x e. 

3.2.1. Strong sticking in pliysisorption. 

In fig|5](left panel), we show the product S x e for surface temperatures of 10 and 125 K, and for 
gas temperatures varying from 100 to 3000 K. This product is the sum of the H2 formation effi- 
ciency times sticking involving physisorbed atoms, Sp/,,., x e/,/,,^, and involving chemisorbed atoms, 
Sf/iem X ^chern - In thesc Calculation s, we consider a sticking coefficient for physisorption, Sphys, such 
as the one in lCuppen et alJ (120101) . and shown in fig|2l For Tj„5,= 10 K, Sxe is very large at low gas 
temperatures and decreases as gas temperature increases. Because of the high sticking probability 
to physisorb, the formation of H2 involves physisorbed atoms and therefore Sxe depends on Sphys- 
At gas temperatures lower than 500 K, H2 is formed through the association of two physisorbed 
atoms (mechanism Ml), while for higher gas temperatures, physisorbed atoms enter chemisorbed 
sites and associate with H atoms in a dimer (mechanism M5). For Trfj„,= 125 K, Sxe is very small 
at low gas temperatures and increases as the gas temperature increases. Because of the high sur- 
face temperature, physisorbed atoms do not stay on the surface and the formation of H2 is insured 
by chemisorbed atoms. The formation of H2 proceeds through direct Eley-Rideal mechanisms in- 
volving a chemisorbed H atom in a dimer (M3), for Tgfli<1000 K, and involving an chemisorbed H 
atom in a monomer for higher Tgas (M2). For the two different surface temperatures considered, the 
mechanisms to form H2 are different since they either involve physisorbed atoms (for Tdusr-^0 K) 
or chemisorbed atoms (for T;/„5,= 125 K). However, the effects of surface dynamics on the forma- 
tion of H2 can be appreciated only by comparing the formation of H2 involving chemisorbed atoms. 
For this purpose we represent in fig|5] (right panel) the different contributions for the formation of 
H2 that involve only chemisorbed atoms; Schem x ^chem- The different mechanisms for the forma- 
tion of H2 that involve only chemisorbed atoms are shown as dashed lines for T;/,„,=125 K and 
solid lines for T^/„.sr=10 K. For Tgas<lOQQ K, H2 forms through Eley-Rideal mechanisms involving 
a chemisorbed H atom in a dimer (M3), while for higher Tgas, H atoms can access chemisorbed 
sites easily even with a barrier, and H2 is insured through Eley-Rideal mechanisms involving a 
chemisorbed H atom in a monomer (M2). The product Schem x Cchem slightly depends on the sur- 
face temperature, but the effect of the surface dynamics cannot be assessed. Indeed, the mechanism 
involving chemisorption depends also on the physisorption since many physisorbed atoms can 
travel on the surface and become chemisorbed. Therefore, as long as some H atoms can become 
physisorbed and subsequently reach the chemisorbed sites, it is not possible to really isolate the 
impact of surface dynamics on the formation of H2. 

3.2.2. Weak sticking in pliysisorption 

In figlH we present the same calculations as previou sly, but we no w consider a very small sticking 



coefficient for physisorption, Sphys, as described in iBuchI ( Il989l) . In this case, Sphys-0.Q2 at low 
gas temperatures and Sphys decreases exponentially as Tgas increases (see fig©. Our results show 
that because atoms do not stick in physisorbed sites, H2 formation is insured by mechanisms in- 
volving only chemisorbed atoms. For Tgfli<1000 K, H2 forms through Eley-Rideal mechanisms 
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involving a chemisorbed H atom in a dimer (M3), while for higher Tg„s, gas phase H atoms can 
enter chemisorbed sites even in the presence of a barrier, and H2 formation is insured through Eley- 
Rideal mechanisms involving a chemisorbed H atom in a monomer (M2). The total product S x £ 
is equal to Schem x ec/,f,„, and the effects of the surface dynamics on the sticking in chemisorption 
can be addressed. If the surface is flat (represented as small dashed lines), then Schem is smaller (see 
fig|2]l and the product is therefore smaller When the surface dynamics are considered, the sticking 
Schem is more efficient, and therefore the product 5 x e is slightly higher. We also note that the 
efi^ect of the surface temperature on the product 5 x e is almost negligible, as shown by the curves 
for T^ust=iO K (solid) and Tj„i,= 125 K (small dashed). 



4. Discussion and conclusions 

In this work, we have studied the influence of the sticking probability on the formation of H2 
using Kinetic Monte Carlo simulations. In our model, mechanisms Ml to M5 (described in the 
introduction) for the H2 formation are included. We have determined the predominant mechanisms 
to form H2 as function of gas and dust temperatures. 

First, we assumed the sticking probability of H atoms to be equal to the transmission coef- 
ficient to cross the barrier against chemisorption. To address the importance of the shape of the 
barrier, we consider a square barrier, and a barrier obtained by DFT calculations. For the latter, 
wavepacket propagation calculations along the coordinate between the H atom and the graphite 
surface (collinear approach) are performed to determine the transmission probabilities of the 
chemisorption barrier. Our results show that if a square barrier is considered, the transmission prob- 
abilities are underestimated by many orders of magnitude at low gas temperatures (Tgaj<2000K). 
The impact of the form of the barrier is visible on the results of the KMC calculations simulating 
the formation of H2. Our results show that the shape of the chemisorption barrier (square or ob- 
tained by DFT) has a strong impact on the H2 formation for Tf/,,,, > 15 K. This impact increases as 
gas temperatures decrease. In conclusion, the formation of H2 calculated with barriers obtained by 
DFT is much more efficient than with square barriers. In these temperatures domain, H2 formation 
involves chemisorbed atoms. For Tj„j(<15 K, H2 is formed by a physisorption mechanism. In this 
case, the form of the chemisorption barrier is not important. 

Second, we investigated the sticking of H atoms in chemisorbed sites taking into account the 
vibrational modes of the surface. For this we used a mixed classical-quantum dynamics method 
where the motion of the H atom i s treated quantum-mechani cally, while the vibrational modes of 



the surface are treated classically (IMorisset & Allouchel2008h . The sticking probability which take 



I 



into account the vibrational modes of the surface can be compared to the transmission probability 
of the barrier The sticking probability is higher by a factor of 5 for T^,,, ~1500K, but is lower 
by a factor 2 (for Tju«=10 K) to 5 (for Tj„s,— 125 K) around Tg„s ~1000K. This calculation show 
that the energy exchange between the atom and the bath of phonon is higher than transmission 
probability of the chemisorption barrier. These sticking probabilities are calculated as function of 
Tg„s, for T(/„s,=10 K and T^„.s,— 125 K. We implement these sticking coefficients Schem in our KMC 
simulations, and show that for T,;„5,=125 K, H2 is formed through a chemisorption mechanism ( 
involving dimers for Tgas< 1000 K and monomers for Tgas> 1000 K ). For Td,„,=10 K, on the other 
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hand, H2 formation involves physisorbed atoms if the sticking in physisorption is high. In this case 
the efficiency to form H2 depends of the sticking probabihties of H atom in a physisorption site. 

The sticking in physisorption influences the formation of H2 at low dust temperatures (<15 K). 
In this stu dy, we considered 2 p ossible sticking coefficients in physisorption Sphys- 1) The sticking 



derived in 



Cuppen et al.l (120101) . which i s valid for me tals and is close to unity for low gas and dust 



temperatures. 2) The sticking derived bv lBuchI (Il989b . whi ch shows that the probability of sticking 



for H atoms on graphite is on the order of a few percent (iLepetit et al. 



2011 . 



Medina & Jackson 



20081). The influence of this Sp/,,., on the product 5 x e is very important at low dust temperatures. 
If the sticking is taken as in 1), then H2 is formed through LH mechanisms with physisorbed atoms, 
and S X e is very high ~ Sphys- On the other hand, if the sticking in physisorption is chosen to be 
as in 2), the formation of H2 involving physisorbed atoms is not efficient, and the formation of H2 
is insured by chemisorbed H atoms. In this case, the product S x e depends on Sci,em x ^cimn- Our 
KMC simulations show that once the vibrational modes of the H-graphite system are taken into 
account in the sticking, the efficiency to form H2 is increased by a factor 2. This behavior is not 
surprising since the sticking probabihty obtained by dynamics calculations is higher than the value 
of the transmission probability of the chemisorption barrier 

In the ISM, the H2 formed are ro-vibrationally excited. The excitation of the molecule is due 
to the energy transfer of the kinetic energy towards the substrate, the translational energy and the 
excitation of the molecule. Th e ro-vibrational excit ation of the molecule allows the determination 
of the mechanisms to form H2 ( ILemaire et alJl201C , 



Islam et alj2007l) . Different mechanisms have 

nf^ — 

been studied theoretically by several groups (IRutigliano et al.l2001 



2003 



2002lMorisset et al. 
2006blBachellerie et al. 



Morisset et al. 



2004: 



1 



Morisset et alj 



Reeetal.12002. 



2005 



Sha & Jackson 



iMartinazzo & Tantardini 



2007n . In our KMC simulations it is not possible to extract the energy of 
the H2 formed. These simulations only allow us to follow the number of molecules formed through 
difl^erent mechanisms (Ml to M5 described in the introduction). In this study, we have determined 
the predomi nant mechanism to fo rm H2 as a function of the dust and gas temperatures. 

Recentlv. lCuppen et alJ (120101) performed similar work on the formation of H2 and the influence 
of the sticking on the product S xe. These authors concentrated on the sticking in physisorbed sites 
(and we used their results in this work), but assumed that the atoms from the gas phase had to 
thermally overcome the barrier against chemisorption to become chemisorbed. In this sense, the 
H atoms coming from the gas phase do not stick in chemisorbed sites at low Tgas {Schem- 0.03, 
0.0001 and lO*^ for Tgas- 500, 200 and 100 K respectively). Therefore, our results are similar to 



Cuppen et al.l (120 lOh for cold dust, when physisorbed atoms are involved to form H2. However, on 
warm dust grains (~ lOOK), chemisorbed atoms are needed to form H2 and the sticki ng of H atoms 



in che misorbed sites sets the formation of H2. In this case, our results diverge from lCuppen et al 



(120 101) for gas temperatures lower than 500 K, since we consider that H atoms from the gas phase 
can tunnel through the important barrier against chemisorption. This process makes the sticking of 
H atoms in chemisorbed sites still efficient at low Tgas- Once these chemisorbed sites are populated, 
H2 can form efficiently through barrier-less routes involving the creation of dimers. Therefore, our 
results show that the formation of H2 is efficient also for intermediate gas and dust temperatures 
(Trf„„>15K and 1000 K>T^„,>100 K). 

Our results suggest that the formation of H2 remains efficient in regions where gas and dust are 
warm. This is true if the time for H atoms to enter chemisorbed sites is shorter than other routes 
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to form H2 (i.e. dense and warm medium). Therefore, the grain surface route proposed in this 
study should be compared to other routes to form H2. Environments whe re grain surface chem 
istry w ould dominate to form H2 includ e photon dominated region s (PDRs, 



Hollenbach & Tielens 



19991) . X-ray dominated regions (XDRs, Meiierink & Spaans 



2004 



2005b . hot cores (ICaselli et al 



1993b 



2009b . Particularly PDRs and XDRs enjoy re- 



and the early universe (iCazaux & SpaansI 
gions of upto 10^^ and 10^"^ cm"^, respectively, where rapid photo-dissociation and chemical re- 
moval of H2 requires it to be fo rmed efficiently on dust grains in order to drive ion-molecule 
chemistry and H2 li ne emission (Habartet al.ll2004b . In PDRs and XDRs gas temperatures can 



be as high as 10 K (iMeiierink et al. 



Hollenbach & McKee 



2007b . values for which the traditional sticking coefficient of 



(11979b decreases rapidly. At early cosmic epochs, for redshifts larger than 
10, the first grain surfaces are expected to be warm, at least 30 K, due to the strong cosmic mi- 
crowave background. Furthermore, the low abundances of metals like carbon an d oxygen cause 
gas te mperatures to be 10~ - lO"*^ K and to be set mostly by H2 and HD cooling (iGlover & Abel 
2008b . Our results indicate that chemisorption effects allow dust grains to act as catalysts for H2 



formation under such hostile primordial conditions. 
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Fig. 1. Transmission coefficients for a square barrier approximation (dotted lines) and for 
rier obtained by DFT calculations (solid lines). Left panel: the different barriers considered, 
panel: transmission coefficients for H and D atoms versus gas temperature. 
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Fig. 2. Left panel: stic king probabili t y in ph ysisorbed sites as function of the ene rgy of t he inc oming 
H atom according to lcuppen et alj J2OI0I) (SI, for Trf„j,= 10 and 125 K) and buchl d 19891) (S2). 
Right panel: trapping probabilities (transmission) and sticking in chemisorbed sites as function of 
the energy of the incoming H atom. The transmission coefficient to overcome a barrier derived from 
DFT calculations is shown as well as the sticking coefficient when the phonons are considered, for 
Trf„j,= 10andl25K. 
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Fig. 3. Left panel: Contour map of the potential energy surface as function of the H position on the 
graphite. The energy difference between two consecutive contour lines is 0.3 eV. A, B, C and D are 
the positions of a barrier between two chemisorption sites. O is the center of the hexagon. H atoms 
can chemisorb and physisorb on top of carbon atoms (A, B, C and D), while they can physisorb 
also in other locations such as a bridge (E). Right panel: stable configuration for the second H 
atom: ortho (closer H atom on the hexagon ), meta (M) and para (opposite H atom on the hexagon) 
configurations and 5 and S (second cycle). dDumont et al.ll2008b 
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Fig. 4. H2 formation efficiencies when squared barriers (61) and DFT barriers (eo) are considered 
to enter in chemisorbed sites. The mechanisms important for the formation of H2 are Ml (2 ph- 
ysisorbed atoms), M2 (direct chemisorption with H monomer) , M3 (direct chemisorption with H 
dimer) and M5 (physisorbed atom entering chemisorbed sites with H dimer) . 
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Fig. 5. H2 formation efficiencies on 10 K and 125 K surfaces as a functio n of gas temperature, 
with the sticking in physisorbed sites S;,/,,., taken as in lCuppen et alj (120101) . Left: total H2 forma- 
tion, governed by mechanisms Ml (2 physisorbed atoms) and M5 (physisorbed atom arriving in 
H dimer) for 10 K dust, and govern by mechanisms M3 (direct chemisorption with H dimer) and 
M2 (direct chemisorption with H monomer). Right; Same as left but for H2 formation involving 
only chemisorbed atoms. At low gas temperatures, H2 is formed through mechanism M3 (direct 
chemisorption with H dimer). For T^aj>1000 K, H2 is formed though mechamism M2 (direct 
chemisorption with H monomer). 
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Fig. 6. Same as figure |5]with the sticking in physisorbed sites Spi,ys taken as in lBuchI (119891) 
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